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Abstract 

p ^ In this paper we compute some of the higher order terms in the asymptotic behavior 

of the two point function P(.4.2(0) < s\,A2(t) < s%), extending the previous work of 
Adler and van Moerbeke [TJ [2] and Widom [25] . We prove that it is possible to represent 
any order asymptotic approximation as a polynomial and integrals of the Painleve II 
function q and its derivative q'. Further, for up to tenth order we give this asymptotic 
approximation as a linear combination of the Tracy- Widom GUE density function fi 
and its derivatives. As a corollary to this, the asymptotic covariance is expressed up 
\ to tenth order in terms of the moments of the Tracy- Widom GUE distribution. 

§1. Introduction 

The Airy 2 process, A2{t), introduced by Prahofer and Spohn [15] in the context of the 
polynuclear growth (PNG) model, is a stationary stochastic process whose joint distributions 
for ti < ■ ■ ■ < t m are given by 

& P (Mh) <s u ..., A 2 (t m ) < s m ) = det (/ - X^2X) L2 ( { i,..., m}X R) (1) 

where = K2(x,y) is a m x m matrix kernel, called the extended Airy kernel, given by 



J °° e t ^ Ai(x + z)Ai(y + z) dz if % > j, 
J ^ e^^i-i-i) Ai(x + z)Ai(y + z) dz if i < j, 



and X = X(x) is the mxm diagonal matrix whose ith diagonal term is the indicator function 
X Si (x) = l x>Si - The one-point function P(yt 2 (£) < s ) is the Tracy- Widom GUE distribution, 

f 2 ( S ) m- ' 

For stochastic growth models and their closely related interacting particle systems, the Airy 2 
process is fundamental since it is expected to describe the limiting process for height fluc- 
tuations belonging to the KPZ Universality Class with droplet (or step) initial conditions 
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(see 0, US] for recent reviews). This has been proved for the PNG model [15], the discrete 
PNG model [13], the boundary of the north polar region of the Aztec diamond [14J and the 
totally asymmetric simple exclusion process (TASEP) [14] . At the the level of the one-point 
function, this universality has been established for the asymmetric simple exclusion process 
(ASEP) [23J and for the KPZ equation [TO HE, El E] • Recent work [16] using replica methods 
have extended this KPZ work to the 2-point function. The Airy 2 process also describes 
the limiting process of the largest eigenvalue in Dyson's Brownian motion model in random 
matrix theory. For further appearances of the Airy 2 process see [6j [7J [10J IT2] . 

We summarize some known properties of the Airy 2 process: 

1. A2(t) has continuous sample paths [TBI fT3] . 

2. A2{t) locally looks like Brownian motion (see Hagg [11] for a precise statement). 

3. The distribution functions ([1]) satisfy nonlinear differential equations [T| |2T 1241. 

4. The covariance cov 2 (t) := cov („4. 2 (t)„4 2 (0)) has the following asymptotic expansions: 

var(F 2 ) -t + 0(t 2 ), t 0+, 
cov 2 (t)={ " C / 1 \ (2) 



n=l 



The small-t expansion of cov 2 was given by Prahofer and Spohn [15] and they also found 
the leading large- 1 term C\ = 0, C 2 = 1. The existence of the higher order terms in the 
large-t expansion of cov 2 was established by Adler and van Moerbeke [TJ [2] and by Widom 
[25] using different methods. In both [TJ [2] and [25J the coefficient C4 was expressed as a 
double integral whose integrand was in terms of the Hastings-McLeod solution of Painleve 
II appearing in the distribution F 2 . One of the main results of this paper is to prove that 
the coefficients C 2n , 2 < n < 5, are expressible in terms of the moments of F 2 . (The odd 
coefficients are all zero.) Precisely, if 



/oo 
s n f 2 (s)ds, f 2 :=K 
-00 



then 



^4 = 2/xi, (3) 
Ce = 2^ + ^, (4) 

13 

C 8 = 2/i 3 + 14^/i! + — , (5) 

126 

C10 = 2// 4 + 24/^3/xi + — // 2 + H6/ii. (6) 

5 

We conjecture that C 2n can be expressed in terms of a polynomial in ^ for i < n — 1. To 
prove these results we follow the program established by Widom [25] and first prove 

n=0 ^ ' 
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as t — > oo. It was previously shown [TJ |2j [25] that for n < 4 each c n could be written as 
polynomials and integrals of the Painleve II function, its derivative, and the variables Si 
and S2- A feature of our analysis is that we show through order t -10 that each C2 n can be 
expressed in terms of fi-, its derivatives, and polynomials in Sk, see (jl8l)-(l20l. ( 1241) . 

Bornemann [U [5] has given a high precision numerical evaluation of cov 2 (t), < t < 100. 
His method involves a numerical evaluation of the Fredholm determinant appearing in ([I]) 
for m = 2 followed by numerical integrations to give cov 2 (t). In Appendix 2 we compare the 
large-t asymptotics with these numerical results. 

In the present paper we begin by showing how to obtain an asymptotic expression for the 
extended Airy kernel following Widom [25]. The large-t asymptotics of the two-point dis- 
tribution is then given in terms of f'2 and its derivatives. This in turn allows the easy 
computation of the large-t expansion of the covariance, which ends the main body of the 
paper. The appendices contain some of the higher order terms, and comparison to high 
precision numerical results. 



The first step in our asymptotic analysis is a large t expression for the extended Airy kernel 
XK2X. In the m = 2 case the XK2X operator has a matrix kernel of the form 



X Sl (x) Ai(x + z)Ai(y + z) dz X Sl (y) X Sl (x) e zt Ai(x + z)Ai(y + z) dz X S2 (y) 



-X S2 (x) e zt Ai(x + z)Ai(y + z) dzXsAv) ^s 2 (x) Ai(x + z)M(y + z) dz X S2 {y) 



We must compute the Fredholm determinant of this operator for large-t. To this end we will 
split the operator into two manageable components so that XK 2 X = K + L: 



2. Asymptotics for xK 2 X 







I 



00 



Xsi(x) 



Ai(x + z)Ai(y + z) dz X si (y) 







K(x,y) : 




■00 







Ai(x + z)M(y + z) dz X S2 {y) 




00 







e z *Ai(x + z)Ai(y + z) dz X S2 {y) 



L(x,y) : 



X S2 (x) / e 



zt A\(x + z)A\{y + z) dz X si (y) 







—00 



3 



The determinant computation is simplified via 

det(J - XK 2 X) = det(J — K — L) — det [(/ - K){I - (I - K)~ X L)] 
= det(J - K) det(J - (/ - K)- l L). 

The determinant of I — K is i 7 2(si)-F2(<S2), and has no dependence on t. So we need only 
look at L to determine the asymptotics. For that determinant, we make an expansion of the 
terms in L(x,y). By repeatedly applying integration by parts, the upper-right corner is 

PCX) 

L12 = X S1 (x) / dz e- zt ki(x + z)M(y + z)X S2 (y) 
Jo 

= E jir E (I) f-Ai« ® Ai<»-">x« + o (t-c«)) 

n=0 fc=0 ^ ' 

and the lower-left corner 

/•O 

£21 = X S2 (x) / e zt Ai(a; + z)Ai(y + z)x si (y) 

J — oo 

= E ^rr- E (I) wu w « Ai, "" ) ^ + ° ■ 

n=0 fc=0 ^ ' 

We take this approximation up to N = 10, which is valid in the trace norm, because this is 
the highest order term we aim to calculate. 

Now we analyze T :— (I — K)~ 1 L. At this point we introduce the notation 

Qn,k = {I ~ X Sk KxiX Sk ) 1 X Sfc Ai ( - \ 

The kernel Km is the classical Airy kernel of random matrix theory, while the purpose of 
the X's are to ensure the kernels are integrated over the appropriate domain. Using this 
notation in the asymptotic expansion of L(x,y), the upper right term in the matrix kernel 
of T is written 

n n / \ 

r » = E ^ E (?)0m « Ai'-*»x, 2 + 0(r<™>) 

n=0 fc=0 ^ ' 

and the lower left term is 

r - = E^S^E ("W®AI<-*>^ + £*-<»«>). 

n=0 fc=0 ^ ' 

The determinant we are now interested in is best analyzed in terms of the following trace 
formula: 

det(J — T) = det(J - T 12 T 21 ) = exp [trlog(7 - T)\ (8) 
= l-trr+i((trT) 2 -trT") -i((tr T) 3 - 3tr T^r T + 2tr T 3 ) H . 
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Here we have introduced T = T\ 2 T 2 \. To evaluate an expression for T we will repeatedly 
use the fact that (/ % g) (e <g> h) = (g, e) L 2f <g) h. By introducing the notation 

u n ,j(s k ) = (Qn,fe,X Sfc Ai (j) ) L 2 = ((5j,fc,X Sfc Ai (n) ) L 2 =u jin (s k ) 

we are able to write 

t= E E ft) n^wflws^n+o 

ni+n 2 <iV-2 fc 1= 0fc 2 =0 ^ ' ^ ' 



(9) 



From (Q above, we know that we will also need the formula for (T12T21) 2 : 

v (-1)"^ ^ ^ ^E ^ / n A / n A / m A / m2 

+mi+m 2 <JV-4 

X M rll -fc 1 ,fc2(s2)Mm 1 -£ 1 ,^(s2)Mn2-fc2A( S l) < 5fei,l ® X Sl Ai (m2 ^ 2) + C (^^T^ ■ ( 10 ) 

We continue multiplying the sums T n in this manner as necessary for the order of our desired 
approximation. 

§3. The Two-Point Distribution 

From OH]) we conclude tr T is order t~ 2 , while fl9]) and (jTOl) tell us (tr T) 2 — tr T 2 is order 
t~ A . Terms can be added as necessary to get any order approximation one wishes for. In 
other words, 

F(A 2 (t) < s 2 ,A 2 (0) <*i) = £ + ° 

for t — > 00. 

It has already been established that c\ = C3 = 0, and c 2 and C4 can be written in terms of the 
Hastings-McLeod solution, which we denote by q{s), to the Painleve II equation (TJ [2J 125] . 
This result can be extended to any order. It is clear from (J9]) and ( ITUi) that all of our 
coefficients c n (si, S2) are polynomials in the variables Ukj multiplied by F 2 (si)F 2 (s 2 ). In our 
notation the first two terms are 

c 2 (si,s 2 ) = F 2 (s 1 )uo fi (s 1 )F2(s 2 )uo,o(s2), (11) 
c 4 (si, s 2 ) = F 2 (si)F 2 (s 2 ) [mi ) o(si)wi, (s 2 ) + ^0,0(^1) (2m 2i o(s 2 ) - Wi,i(s 2 ))] (12) 
+ reversed, 

where "reversed" denotes the interchange of si and S2 in the previous terms. The formula 
for F 2 in terms of q, 

r r°° 

F 2 (s) = exp — / (x — s)q(x) 2 dx , 

is well known so we need only focus on the Ufc,/s. There are three facts that will enable us 
to show each Ukj can be expressed in terms of q and q' and integrals of q and q'\ 
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1. go : = (I — Km) 1 X a Ai(s) is a solution to Painleve II (q = q), 

2 - u 'k,j = -QkQj, where qj := (I - K Ai )~ 1 X s Ai (j) (s), 

3. q n = (n- 2)g„_ 3 + sq n ^ 2 - m«-2,i q + «n-2,o Qi for n > 3. 

The first statement is established explicitly in [20j 122] , where a template is also developed for 
the proof of the second statement. We begin by proving the second statement. Henceforth 
we will use the notation 

Q k (x;s) := (I-Km^XsAiW. 

Proof of Statement 2. We proceed by first differentiating u k j- 
duk j d 



Qk(x; s)Ai^\x) dx 



q k (s) Ai® (a) + ^Ai^(x) dx. (13) 

J s OS 



ds ds 
We now need a formula for d s Qk- 

^ = ^(1- KmT'XsA^ = (I- K M y^(I - K M )^ Xs A^. 

OS OS OS 

This necessitates a formula d s Kx\. Direct computation shows that this has kernel 

— = -K M (x, s)5(y - s). 

By introducing the operator R = (I — K^)~ x K^x we can write 

dQ k 
ds 

Returning to (ITBl . we now have 



(x; s) = -R(x,s)q k (s). 



^ = -q k (s)A^(s) - [°° R(s,x)q k (s)Ai«\x) dx 



(5(s -x) + R(s, x))Ai {k) (x) dx 



= - qj (s)(I - K^XsAi^is) = - qj (s)q k (s). 

Note that we pass from the the second line to the third line because the kernel of (/ — Km) 1 
is 8(x — y) + R(x,y). □ 

Proof of Statement 3. First, we develop an explicit formula for q n . We start with 

Q n (x-s) = {I-K M )- l X s Ai^ 

and the formula 

Ai (n) (x) = (n- 2)Ai (n - 3) (x) + xA\ {n ~ 2 \x) 
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for n > 3. This gives 

Q n {x; s) = (n- 2)Q n „ 3 (x; s) + P n - 2 (x; s) 
where P n (x; s) := (/ — Kp Li )^ 1 xAi^ n \x)X s - Now we apply the commutator relations 

[M, (I - Km)- 1 ^ = xQ n (x; s) - P n (x; s) 

and 

[M, {I - Km)- 1 } = Q {x; s)Q 1 (y; s) - Q 1 {x; s)Q {y; s) 

to arrive at 

P n (x; s) = xQ n (x; s) - u n ,i(s)Q {x; s) + u n , (s)Qi(x; s). 

Taking x = s we get the recursive relation 

q n = {n- 2)g„_ 3 + sq n - 2 - «n-2,l? + «n-2,ogi- 

for n > 3. From [2H E21 EH] we know that qi = q' + «oo 9 and q 2 = sq± — «u q + Moo <Zi, so we 
have formulas for all q n . □ 

From here, we can complete the proof inductively. If Ukj is expressible as a polynomial and 
integrals of the q, q' and s for j < k then so is qk+i due to the identity above. The differential 
equations for each Ukj give us Uk+ij = qk+i(x)qj(x)dx for j < k + 1. This completes the 
proof, so long as enough base cases are satisfied. 

While it is interesting that we can write these Ujk formulas as a polynomial and integrals in 
terms of s, q, and q', there are explicit polynomial formulas in the aforementioned terms for 
each Uj^k with j + k < 8. For example, 



Uq, = 


{ /\2 2 4 

(q) - sq - q , 




(14) 


Ulfl = 


1 2 1 2 

2*0,0-29. 




(15) 


Ul,l = 




l 2 ' 

I M o,o - 3 99 > 


(16) 


«2,0 = 


1 2 1 

2 U i,i + 2 




(17) 



We conjecture this is true for all Uj k- The conclusion is the two-point distribution has an 
asymptotic formula up to order t -10 in terms of q. 

The formulas for q n and u'j k also enable us to show that F 2 Uj^ is multilinear in f 2 = F 2 , 
its derivatives, and s for j and k needed through ciq. A table of these formulas is included 
in Appendix 1. The verification of each equation simply requires the differentiation of both 
sides and looking at s — >■ oo to verify the constants of integration are zero. 
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Using this table, (J9J), and ffTOj) . a direct computation yields the formulas 

c 2 (s u s 2 ) = / 2 (si)/ 2 (s 2 ), (18) 
c 4 ( Sl)S2 ) = (si + s 2 )/ 2 ( Sl )/ 2 ( S2 ) + ^( Sl )^(s 2 ), (19) 

c 6 (si, s 2 ) = i (3si + s 2 )(3s 2 + si)/ 2 (si)/ 2 (s 2 ) + 3 (/ 2 (si)/ 2 (s 2 ) + ^(si)/^)) 

+ s 2 )^( Sl )/i( a3 ) + i tt{si)tt{s 2 ). (20) 
o 

The odd terms are not listed here because they are zero. The coefficient formula for cs is 
given below in ( 1241) . 

From a logical standpoint the work we have done is sufficient. However, it is instructional 
to see an alternate derivation of (fIB]) - (JT9~]) from (fTT]) - (fl2"]) . First note that we have already 
established wo,o(s) = / q(x) 2 dx. Using this and the integral formula for F 2 , it is easy to 
check that / 2 = F' 2 = F 2 m 0j o- This gives us (JT51) . For the next term, we use the identities 
( 1T7|) . (TT5T) . and / 2 = F 2 u fl in the formula (jl2"l) to get 



c 4 (si, s 2 ) = - F 2 (s 1 )(u fi(s 1 ) 2 - g(si) 2 )F 2 (s 2 )(M ,o(s 2 ) 2 - q{s 2 ) 2 ) + (si + s 2 )/ 2 (si)/ 2 (s 2 ). 
Finally, observe that = [-^2^0,0]' — -^2^00 — -^2<? 2 so we arrive at (jl~9l) . 

§4. The Covariance 

For the covariance of the Airy 2 process we denote the iV-th order asymptotic approximation 
by 

N c 

cov 2 , N (t) = J2~^- 

n=l 

It was established by Adler and van Moerbeke [H [2] and Widom [25J that the coefficients up 
to iV = 4 are d = 0, C 2 = 1, C 3 = 0, and 



C 4 = / / c 4 (w,f) dudv. 
J Ju 2 

Through an application of Fubini's theorem and (|T9|) . we get 



1 ' - 



C A = 2 / uf 2 {u)du / f 2 (v)dv + -[ I f 2 {u)du) =2^. (21) 



Following this same procedure for Cq we have 

C G = / c 6 (s 1 ,s 2 )ds 1 ds 2 
Jr 2 



/ ~ (3si + s 2 )(3s 2 + si)/ 2 (si)/ 2 (s 2 ) dsi ds 2 
Jr 2 6 

2/i 2 + y/i? (22) 



since all the other integrals appearing in fj20|) integrate to zero. The analogous formulas for 
C 8 and Cio appear in (jSJ) and (Q with some details of their computation in Appendix 1. 
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Appendix 1: Identities 



In this appendix we give the formulas necessary to attain f|T8|) - fl20|) . The table below gives 
Uj,kF2 in terms of f% and its derivatives. 



j, 


K 


h 


-f ' 
/2 


I'll 


f(3) 
J2 


J2 


/(5) 
J2 


f (6) 
J2 


A7) 

J2 


As) 
h 


0. 





1 


























1. 








1 

2 





{-\ 



















1. 


1 


g 

~3 





1 
3 




















2. 





g 

3 





1 

6 




















2. 


1 


1 

_ 4 








1 

8 

















3. 





7 

1~2 


g 

3 





1 

24 

















2. 


2 


s 2 

5 


3 

— To 








1 

20 














3. 


1 


s 2 

~y 


2 
15 


6 





1 

30 














4. 





s 2 

y 


17 
60 


6 





i 

T20 














3. 


2 


2s 
~9 


s 2 
18 


1 

~V2 


Ts 





1 

72 











4. 


1 


13s 

~~ Ts 


s 2 
_ 18 


11 

24 


g 

9 





1 

T44 











5. 





101s 
90 


23s 2 
90 


59 
120 


Ts 





1 

720 











3. 


3 


s 3 34 

— y ~"~ 63 


17.s 

yr 


s 2 

y 


i 

_ 18 


36 





1 

252 








4. 


2 


s 3 11 

y — 42 


13s 
_ Y4" 





1 

8 


2T 





1 
336 








5. 


1 


s 3 74 

— y ~~ To5 


47s 
_ 420 


s 2 

To 


7 

20 


2T 





1 
840 








6. 





s 3 i 1151 
~ + 630 


733s 
420 


7s 2 
45 


71 
360 


g 

72 





1 

5040 








4. 


3 


s 2 
4 


127 

288 


5 s 
Ts 


s 2 
18 


5 
288 


g 

72 





1 

1152 





5. 


2 


43s 2 
60 


s 3 123 
15 _ 160 


3s 
_ 20 


s 2 
30 


61 

480 


60 





1 

1920 





6. 


1 


73s 2 
60 


s 3 _|_ 271 
15 1440 


29s 
36 


17s 2 
180 


221 
1440 


s 

90 





1 

5760 





7. 





691s 2 


22s 3 _|_ 4873 
105 1440 


394s 


lis 2 


83 


s 





1 





420 


315 


180 


1440 


360 


40320 


4. 


4 


s 4 118s 
9 81 


17s 2 
54 


s 3 | 119 
81 1 144 


31s 
108 


s 2 
27 


1 

144 


s 
216 





1 

5184 


5. 


3 


s 4 _|_ 667s 
9 ^ 810 


197s 2 
540 


32s 3 4 
405 45 


119s 
1080 


29s 2 
1080 


11 
360 


lis 
2160 





1 

6480 


6. 


2 


s 4 , 322s 
9 i " 405 


103s 2 


s 3 35 


113s 


37s 2 


11 


s 





1 


540 


162 72 


540 


1080 


180 


216 


12960 


7. 


1 


s 4 21167s 


1999s 2 


37s 3 , 115 
567 ' 63 


6107s 


47s 2 


17 


s 





1 


9 5670 


3780 


7560 


1080 


360 


432 


45360 


8. 





s 4 | 19912s 


5297s 2 


409s 3 | 28319 


4273s 


19s 2 


19 


s 





1 




9 ~ l ~ 2835 


1890 


2835 1 10080 


7560 


1080 


1440 


2160 


362880 



10 



Thus, for example, u 2j ±F 2 = — I/2 + §/•} ■ To verify these identities, first divide both sides 
by F 2 , differentiate both sides, then compare the behavior at infinity to check that the 
constant of integration is the same. To illustrate this, we check the identity for Ui^F 2 . After 
differentiation we must verify that 



-?i9o 



F 2 



which is equivalent to 



-(<?' + u 0fi q)qF 2 = - (-2qq'F 2 - 2q 2 f 2 + 1*0,0/2) ~ \ f 2 u o,o- 

Expanding each side, we see they are both —qq'F 2 — q 2 f 2 . The remainder of the table is 
proved in this manner. 

For n = 2,4, 6, the identities presented in this table are sufficient to write c n in terms of f 2 
and its derivatives. At c§ we have a potential problem. There is a nonzero eighth order term 
appearing in | [(tr T) 2 — tr T 2 ]: 



f 8 



[«i,o(si) 2 - «o,oOi)«i,i(si)] [mi j0 (s 2 ) 2 - Wo,o(sa)«i,i(s2)] 



This term must be multiplied F 2 (si)F 2 (s 2 ), and the table does not provide a way to reduce 
this to an /2-type expression. However, a direct computation will verify the formula 



[u 2 j0 - U ,0«l,l] F 2 



-lf2 + \sf 2 - 1/. 



(iii) 



6" 3 " z 12 

Using this formula, the table of identities, and © and ( ITU]) we arrive at the expression^ 



(23) 



cs(si, s 2 ) 



+ s? + 7S1S 2 . + 7s 2 s 2 + s 3 2 



149 

6 

34 

15si + — S 2 



3 2 13 :.) ■■■) 

— Si T" Si So \ ~ So 

2 1 3 2 2 



3 



f2(si)f 2 (s 2 ) + 



34 

15s 2 + —si 



Performing the same sort of integration as in (T2"Tj) . (|22|) results in (jSJ). 



c(iii) 1 



(24) 



We can find t~ 10 order terms in trace formulas of the form tr T J (tr T) k for j + k < 5 in the 
expansion (jSJ). Remarkably, when the computation of tenth order terms in (jSJ) is carried out 
it is found that all terms cancel except those contained in 

-tr T+i [(tr T) 2 -tr T 2 ] . 



1 We do not explicitly provide the formula for cio, as the result is rather cumbersome. 
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The terms needed from — tr T can be dealt with using the table, as before. The same is not 
true of |[(tr T) 2 - tr T 2 ]. There is a great amount of cancellation within this term, so much 
so that by introducing the notation 

A{s) = u i,o( s ) 2 ~ Mo,o(s)ui,i(s) 

B(s) = -2w 2 ,oO) 2 + wi,i(s)m 2i o(s) - ui t o(s)u 2 ,i(s)2uo fi (s)u 2 ,2(s) 
+3u li0 (s)u 3i0 (s) - 3u ,o(s)«3,i(s) 

= «l,o(s)«2,o(s) - Mo,o(s)n 2 ,l( s ) 

we find the tenth order term in |[(tr T) 2 — tr T 2 ] is given by 

2A(s 1 )F 2 (s 1 )B(s 2 )F 2 (s 2 ) + 2A(s 2 )F 2 (s 2 )B(s 1 )F 2 (s 1 ) + 2C( Sl )F 2 ( Sl )C(s 2 )F 2 (s 2 ). 

Analogous to the above table and (|23|) . we may verify the following identities in the same 
fashion: 

A(s)F 2 (s) = -i/ 2 ( a ) + l a ^( a )_±/W( a ) 

B(s)F 2 (s) = -lsf 2 (s) + ls 2 f>(s)-±f!?\s) 

C(s)F 2 (s) = ^m + lsf^-^is). 

Using these formulas to integrate terms in |[(tr T) 2 — tr T 2 ], and the table of identities to 
do the same for terms in tr T, we find for C\q. 

Appendix 2: Numerical Comparisons 

Folkmar Bornemann has generously provided us with a table of values for the Airy 2 co- 
variance obtained by directly computing appropriate Fredholm determinants. See (4J \S\ for 
details on the numerical methods used. We denote these values cove(t) and consider the 
information for 5 < t < 25. Using the following high precision computations courtesy of 
Prahofer, 

Aii = -1.771086 807411601626... 
/i 2 = 3.949 943 272 220 377 513... 
^3 = -9.711844 753 027647354... 
/i 4 = 26.025 435 426 839 994 565... 

together with Q-©, we obtain the expressions 

C 4 = 2/ix = -3.542173 614 823 203 252... 

C 6 = 2/i 2 + —fil = 18.355 714 809 065 487... 
3 

13 

C 8 = 2/i 3 + 14/x 2 /ii + y = -110.863 383 378 407421... 

126 

Cio = 2/x 4 + 24// 3 /ii + + 116/xi = 652.588 990 733 866 004 •• • . 

5 

12 



The following table collects values of covs(t), the approximations cov2 5 „(t), and the error 
between the two values measured by cove(i) — cov 2jn (f). 



t 


cov B (t) 


COV 2 ,6(£) 


Error 




COV 2 ,s(i) 


Error 




COV240W 


Error 




5 


.03527955721 


.03550728796 


-2 x 10 


-4 


.03522347770 


6 x 10" 


-5 


.03529030281 


-1 x 10" 


-5 


10 


.00966309240 


.00966413835 


-1 x 10 


-6 


.00966302972 


6 x 10" 


-8 


.00966309498 


-3 x 10" 


-9 


15 


.004376044913 


.00437608706 


-4 x 10 


-8 


.00437604380 


1 x 10" 


-9 


.00437604493 


-2 x 10" 


11 


20 


.002478143955 


.00247814822 


-4 x 10 


-9 


.00247814389 


6 x 10~ 


11 


.00247814396 


-1 x 10" 


12 


25 


.001591006500 


.00159100722 


-7 x 10" 


10 


.00159100649 


6 x 10" 


12 


.00159100650 


-8 x 10" 
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